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(57) ABSTRACT 
CTEeiinyeln^^ 

^c^fficientsT to^tT^ule~dby"a tim edomam'eq uaH^tionxircuiC 
i£]Tre^ 

multitone(DMT)~technologyrThe method and apparatus is 
particularly efficient in that the order of the channel response 
equation after channel shortening is less than or equal to the 
allotted cyclic prefix length utilized in the applicable com- 
munication protocol being used. Further, the invention 
includes an efficient method and apparatus for determining 
symbol synchronization for the time domain equalization. 
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METHOD AND APPARATUS FOR TIME- 
DOMAIN EQUALIZATION IN DISCRETE 
MULTTTONE TRANSCEIVERS 

FIELD OF INVENTION 

Hie invention pertains to transceivers and modems. More 
particularly, the invention pertains to asymmetric digital 
subscriber line (ADSL) modems used to achieve very high 
speed data communication via telephone networks. 

BACKGROUND OF THE INVENTION 

There is an ever present desire to maximize the speed of 
digital communications via networks, and particularly tele- 
communication networks, such as public telephone systems. 
Accordingly, telecommunication network providers now 
offer to their customers many options for coupling to the 
telephone network in addition to the standard analog based 
connection commonly referred to as POTS (Plain Old Tele- 
phone System). Some of the options that are widely avail- 
able are integrated services digital network (ISDN), T-l 
lines, E-l lines, digital subscriber lines (DSL) and asym- 
metric digital subscriber lines (ADSL). ADSL/s can provide 
very high data speeds such as on the order of several 
megabits per second, over a standard twisted wire pair. 
Unlike the traditional data modems used for analog com- 
munication with a telephone central office via a twisted wire 
pair, ADSL requires modems both at the subscriber end and 
at the telephone company Central Office end. Current ADSL 
systems employ discrete multitone (DMT) technology to 
implement high bandwidth communications, such as for 
digital TV broadcast, on demand video, high speed video 
based internet access, work at home digital file transfer, 
teleconferencing, home shopping, and information services 
over existing twisted wire pair telephone lines. 

The international telecommunications union (ITU) has 
promulgated a standard for ADSL that is commonly termed 
G.lite and which is set forth in ITU-T specification G.992.2, 
incorporated herein by reference. Another standard promul- 
gated by ANSI is commonly termed Heavy ADSL and is set 
forth in ANSI specification T1.413, issue 2, also incorpo- 
rated herein by reference. G.lite has 256 samples with 128 
tones (32 tones for upstream communications), each tone 
having a real and imaginary portion. G.lite uses a cyclic 
prefix of length 16 samples. Heavy ADSL has 512 samples 
with 256 tones (32 tones for upstream communications), 
each tone having a real and imaginary portion. Heavy ADSL 
utilizes a cyclic prefix of length 32 samples. 

FIG. 1A is a block diagram of the basic ADSL modem 
functions in accordance with ITU-T G.992.2 with 256 
samples per symbol. The upper half of the diagram repre- 
sents functions in the transmit direction while the lower half 
represents functions in the receive direction. 

In the transmit direction, digital data is generated in block 
102. That data is processed through a scrambler 104 and then 
through a forward error correction (FEC) encoder 106 which 
adds syndrome bytes to the data that will be used for error 
correction by the receiver at the receiving terminal. Next, as 
shown in block 108, the transmit data is encoded using 
quadrature amplitude modulation (QAM). The data is then 
converted from the frequency domain to the time domain via 
inverse fast fourier transform (IFFT) 110. 

A 1:4 interpolator 112 interpolates the output of IFFT 
block 110 to produce 256 samples from the 64 samples 
output from block 110. Acyclic prefix is added to each frame 
in block 114. The cyclic prefix comprises data added 
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between the symbols to avoid inter symbol interference 
(ISI). The cyclic prefix is of a standardized length, e.g., 16 
samples for ITU-T G922.2 or 32 samples for ANSI-T1.413, 
issue 2, and is removed at the receiver to recreate the original 
5 transmitted data. The data is then forwarded to a coder/ 
decoder (CODEC) 116. The CODEC encodes the data for 
transmission over the twisted wire pair to the receiving 
device. 

<SlSefiiy?^^r©>n of the transceiver 100, the received 
10 signal is passed from the twisted wire pair through the 
CODEC 116 where it is decoded. I t|isltH^p.asse:d f 4to;a£ Um_e3 
^dpgaii^^ 

^sp^n se>Thenrirri2Q, the cyclic prefix is removed. Next, 
an echo canceller (EC) 134 creates an echo cancellation 

15 signal based on the transmit signal which is subtracted by 
subtracter 121 from the receive signal in order to cancel any 
echo of the transmit signal that might return over the twisted 
wire pair and interfere with the receive signal. The echo 
compensated signal is converted back to the frequency 

20 domain by fast fourier transform (FFT) in 122. T!rielpinfl24? 
^e^uene.y^domai^ 
pjfeateflfor^tr^ 

en^et(jISI)) The receive signal is then processed through a 
quadrature amplitude modulation decoder (QAM) 126 to 

25 decode the tone signal into digital data. That is followed by 
forward error correction (FEC) 128 which uses the syn- 
drome bits that are added by the transmit path FEC encoder 
106 to perform forward error correction. Finally, the data is 
descrambled to extract the true data signal in 130 and then 

30 forwarded to the receiver 132. 

FIG. IB is a block diagram of the basic modem functions 
in accordance with ANSI-T1.413, issue 2. Note that the 
functions are similar to those for ITU-T G.992.2. Notable 
differences include that the receive path FFT 122a is a 512 

35 (rather than 256) point FFT and the transmit path interpo- 
lator 112a is a 1:8 interpolator (rather than a 1:4 
interpolator). Both differences are a result of the different 
number of samples between the two protocols. 

FIG. 2 is a block diagram which helps illustrate time 
domain equalization in DMT in more detail. In particular, 
block 202 encompasses functions performed at the trans- 
mitter. Block 207 encompasses functions that occur in the 
data channel and block 212 encompasses functions per- 
formed at the receiver. Block 206 encompasses those func- 
tions which result in channel shortening. Thus, it can be seen 
that, in the transmitter 202, a DMT signal, x(n), is generated 
f or tran smission via the channel. /H^pnl&®lf207JKai^ 
>cffa^nel>ffi 
A^elirclutgiy^ 
the^cony^p^ 
^^sppnse^h(n), i.e ; 

M»M?)V&P Eq. (1) 

flflit^tim^d blaaTfe 
55 21&Jeonwly^ 
l^(rl)^3foe~j o^p^ 
^oSffi:circuifell8^ therefore^is : 

Eq. (2) 

60 ^^rEQxif^iQ18f shortens the criMnel r^the"lengthrof ^ 
thT^yclic'prefix^e 

HcavyjADSL. Then the cyclic prefix is removed from the 
signal y'(n), as shown in block 120 in FIG, l.^^c^hj(n)' 

65 ^annel^theM ' ^~ 

eniggg t (^^ 
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tbL^.OLequalJojE^le^^ FIG. 4 is a flowchart illustrating symbol synchronization 

itHen eacrTy'(h]i^ffl"^ve.no^epe ndcncy on _anyJransmitted in accordance with the present invention. 
<sampleA_bef ore-x(n-L) . 

Afiej^thT^liOre^ DETAILED DESCRIPTION OF THE 

(^athouF inter symbol interference) shou ld„remajnl 5 INVENTION 

Howevej:,jLM^I^the^^ 

dent onjthejransmit samphng-sigDal-that-appears-before 1 r The goal of the present invention is to design a TEQ ^ 

.ffgL^d ^ shortening filter that produees_a_new^^ 

pr^ent~~inter-sym^^ [respo nse h'(n ) forwhi ch MgL . 

performance_of „the_transceiver._Hence,-it^is-desirable-to 10 The impulse response h(n) of the channel has the transfer 

^designWefficient TEQ shortening fiKer;j^a(n Vthat make s function : 
MgU^nd^fdFwfaich J he enef ly^ inZthQhamellimMlse 

(response r^(n)Jsjis_^^ ~ ^--bo"+-biz- l -+ :t>fz' 2 + a u ~7 Et i- ( 3 ) 

^^ssible for^<LT : ^ - - ^ 

SUMMARY OF THE INVENTION 



15 



g^ny^nto^-a^ne^ where L ^ the length of me cydic prefix> K is me ]ength of 

^hejyeffi cientsno-bc^u^ the TC q filter having impulse response (l+a(n)) or A(z). In 

cir^ Jnj j^geiy^jn_a commuriica^o^ystem Jising a accordance with equation 3, the transfer function of the 

< disj^te^ultit_on^ 20 shortened channel is: 
^t^^particinaTl^fficien^ 

^gon^_e^uauQn_aftel3ha //(z)^A(r)//(2)=fc 0 +6 1 2- 1 +i 2r - 2 + . . . +fc (i _ 1 >r< J -- 1 > Eq. (4) 
t equ^~ to~the~allottedn cycfo 

^omrmm^don protorolrFurther, the invention includes an Thus, the length, M, of the shortened impulse response shall 

efficient method~and apparatus for determining symbol 25 be set equal to the length of the cyclic prefix, L> in order to 

synchronization for the time domain equalization. eliminate ISI between two consecutive symbols by remov- 

More particularly, the invention is a method and apparatus tne cvcuc P refix to improve the overall performance of 

for generating coefficients for performing time domain me transceiver, 

equalization on a digital multitone signal involving (1) Let 

synchronizing the received DMT signal to a corresponding 30 

signal that was transmitted from said transmitter prior to ■ ■ ■ Vir fl i ~ a 2 • ■ ■ - a (K-l)] ^ -[&-<^f Eq. (5) 
inclusion of any channel response, (2) <detennining-time 23 
v do^ainreq uaIization~coefficients "vi a^ 

(ff^^ -in _r ^ ( n -i r x_ 35 W(«Mx(k>(«-1Mk-2) ■ ■ • 'HW^NM""!) ■ • • 

. ^ ■ l)f Eq. (6) 



or 



where T denotes the transpose of a matrix. From equations 

^^^^^r^JJ&rf^^ii)- (3) through (6), the estimated signal, y(n), at the receiver is: 



where 



?{n)=I T W{n) . Eq. (7) 

and the error is: 



$^^{y(n)y(n+k)}=Ry(-ky «^)^K»H(") Eq. (8) 

45 

r f i From the well known LMS (Least Mean Square) error 

^ ygffl?.) w t"Jr- 1~ J'o algorithm, the optimal solution 1^ is 



and 50 
T represents a transpose function; and 



RI ep rr Eq. (9) 

where 



(3) time domain equalizing the received signal using said R-E{w(n)w r (n)} Eq (10) 
determined coefficients. 

BRIEF DESCRIPTION OF THE DRAWINGS *5 ^{y("M»)} Eq. (n) 

FIG. lAis a block diagram illustrating the functions of an E{ } in the equations above represents a generalized aver- 

ADSL modem transceiver in accordance with ITU-T G ag^g function or an ensemble average. \ opl must be gener- 

992.2. ated before a value for a can be found. 

FIG. IB is a block diagram illustrating the functions of an 60 Let 
ADSL modem transceiver in accordance with ANS1-T1. 413, 

issue 2. ^)^W»M"+fr)}^xrK) Eq. (12) 

FIG. 2 is a block diagram illustrating in more detail R x m=E{ x ( n )y(n + k)}=R yx (-it) Eq . (13) 
channel shortening in discrete multitone communications. 

FIG. 3 is a block diagram illustrating the functions of an 65 * y #>E{y(«M"+fr)}=*x>>(-*) Eq. (14) 
ADSL modem transceiver for digital multitone communi- 
cations in accordance with the present invention. With equations (9) to (14), we have 
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Equation (15) can be rewritten as 



20 



Eq.(2S) 







r t> opt 









where 



iWO) 

fed) 



*W(2) 



Eq. (16) 



25 



/?jy(2) 



« v (-2) 

M- 1 ) V-2) 



K v (L-3) 



*^(L-2) 

. /?vy(A:-3) 
■ *j»(0) 

. /? XJ (-(^-2)) 



Eq. (17) 



Eq. (18) 



Eq. (19) 



Equation (25) can be solved in several ways. In order to 
reduce computation time and memory usage, some condi- 
tions will be considered and some matrices can be pre- 
calculated and stored in a table. In equation (25), is an 
autocorrelation of the training signal x(n) that is pre-defined 
in ITU-T G. 992.2 and is a positive definite symmetric 
matrix. Hence, the matrix (R^R^R^-R^) also is a 
positive definite matrix. The implementation for equation 
(25) can be summarized as follows: 

(A) Since R^" 1 is a positive definite matrix, it can be 
rewritten as R ;c;c ~ 1 =V r D £rIV V, where V r is a lower triangular 
matrix whose main diagonal elements are all l's, and D t -„ v is 
a diagonal matrix. These matrices can be pre-calculated and 
stored in a table. Only L(L-l)/2 locations are needed to store 
the V matrix and L locations are needed to store the D /m , 
matrix. Hence, equation (25) becomes 



40 



and 



Eq. (20) 45 



Eq. (21) 



55 



In the equations above, R^ and R yy are positive definite 50 
symmetric matrices with the dimensions LxL and (K-l)x 
(K-l), respectively. R^, r 2 and r 2 are Lx(K-l), Lxl and 
(K-l)xl matrices, respectively. We are interested in finding 
the coefficients of A(z), that is a opf . Hence, from equation 
(16), we have 



R^^-R^^r, Eq. (22) 

RjKp-Ryy^opfr z . Eq. (23) 

From equation (22), we can solve b opt as 

K p rR^{r^R xy a 0 ^ Eq. (24) 

By substituting equation (24) into equation (23), we have 

{R Xy ^j^" x R X y-Ry^ < ^T2-QL x *R. X] ^ x rd 

or 



(B) The matrix VR^ which appears 3 times in equation 
(26) can be computed once and used in all three instances. 
The matrix F^VR^D^ VR^-R^,) in equation (26) is a 
(K-l)x(K-l) positive definite matrix. 

(Q By employing the Cholesky decomposition method 
(L. R. Rabiner and R. W/ Schafer, "Digital Processing of 
Speech Signals", Prentice-Hall, Inc., Englewood Cliffs, N.J., 
1978), equation (26) can be solved for a 0 ^ r 

Equation (3) is a general formula that covers cases of 
causality than can occur in the real world. Depending on the 
behavior of a particular design, the transfer function may 
alternately be written as: 



A(z) do +a 2 z~ 2 + ... +o t -i ) r ( * -1 * 



Eq. (3a) 



60 



65 



and the error e(n) that is to be minimized may be written as: 

*(«M«H(«) Eq.(&) 

where x is the estimated signal at the receiver. The same 
concepts described above to arrive at Eq. (25) based on Eq. 
(3) can be used with respect to Eq. (3a) to derive a similar 
solution. 
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FIG. 3 is a block diagram illustrating the operation of a 
modem transceiver for digital multitone communications in 
accordance with the present invention. It is essentially 
identical to FIG. 1, except for the new TEQ filter 302 
implementing the operation described hereinabove. s 

To effectively implement the new computation of TEQ 
coefficients in accordance with equation (25) or (26), y(n) 
from the CODEC should be well synchronized to the 
transmitted signal x(n). If the received signal y(n) is not well 
synchronized to the transmitted signal x(n), then R^(k) in 10 
equation (15) may not generate the correct value and the 
TEQ coefficients obtained would not be an optimal solution. 
Accordingly, it is important to accurately establish symbol 
synchronization to achieve the best performance of the 
above algorithm. 15 

A technique for symbol synchronization will now be 
discussed. Let us consider handshaking between two DSL 
modems. The Central Office (CO) modem sends the 
C-REVERB2 signal (see ITU-T G.992.2 specification for 
definitions of variables) to the subscriber's modem. The 20 
subscriber's modem uses the C-REVERB2 signal to perform 
symbol synchronization and compute the TEQ coefficients 
simultaneously. In order to find the symbol synchronization, 
a cross correlation technique is employed as 

25 

N " N Eq. (27) 

n 4- m ~ k) 



* V* x{n)x{n 4- m - k) 

N 



30 



where N is the number of samples in one DMT symbol, 
0^^) is the cross correlation of x(n) and y(n+m), and 
0xJLm) is the autocorrelation of x(n) and x(n+m). Since FFT 
(Fast Fourier Transform) is an indispensable module in a 40 
DMT DSL modem, we can save memory and computational 
overhead by calculating 0^(01), O^m^N, in the frequency 
domain (i.e., FFT). By converting both sides of equation 27 
by FFT, we derive: 

<D^(o)^/T(^H)-if(co)^(a)) Eq.(28) ^ 

where H(co) is the FFr of the channel response h(k) Let x(n) 
be the C-REVERB2 signal and P(co)=FFT(x(n)), Then equa- 
tion (28) can be written as: 



O v ((o)=i/(o>)<I> JCK (a))=HC6i)P(ai)P*(a>)=y(a))P'((o) 



& r (m)=rFFT(<I> xy {u)) 



50 



Eq. (29) 



where P*(w) is the conjugate of P(co) and Y(co)=FFT(y(m)) 
in equation (29). From equation (29), the cross correlation of 
0 xy (m) for O^m^N can be easily found by employing an 5S 
IFFT (Inverse Fast Fourier Transform) module as 



Eq. (30) 



0^(m) for O^m^N in connection with a subscriber's DSL 
modem also can be interpreted as the impulse response of 60 
the whole channel, including (1) the transmitter filter of the 
CO DSL modem, (2) the characteristics of the two wires 
between the CO modem and the subscriber's modem, and 
(3) the receiving filter in the subscriber's modem. 

There are many factors that determine which sample 65 
between 0 and N should be selected as the start point. 
Therefore, the cross correlation R^ and the auto -correlation 



R yy in equation (15) should be computed. Two major factors 
in determining the symbol synchronization are the length of 
the TEQ filter and the characteristics of 0^,(m). The length 
of the TEQ filter determines how many peaks and valleys of 
0^(01) will be required to get the desired performance of the 
TEQ coefficients. The characteristics of 0^(m) determine 
the spacing between the peaks and valleys. 

A flowchart to implement the symbol synchronization in 
accordance with the description above is shown in FIG. 4. 

In step 402, the impulse response of the channel is 
determined by utilizing FFT techniques or cross correlation 
techniques. In step 404, the impulse response ^(m) for 
O^m^N passes through an energy detector with a window 
of the size of the channel prefix length, termed CPLEN, in 
order to determine the maximal energy and the range of the 
channel signal. In step 406, the noise floor of the channel is 
determined from the entire channel impulse response. The 
average noise floor can be calculated easily by avoiding the 
signal energy around the maximal energy area. 

In step 408, the start point of the channel signal is 
determined by using the noise floor condition. Then in step 
410, starting at the start point determined in step 408, the 
values of the first C peaks and valleys that are totally within 
the channel impulse response are determined, where C is a 
pre-determined constant based on experiment. For example, 
in the implementation of ITU-T G.992.2 for ADSL modems, 
CPLEN is 16 and the length of the TEQ is set to 6. With 
these conditions, the number of peaks and valleys, C, can be 
set to 6. 

In step 412, equation (25) is used to solve for the TEQ 
coefficients starting at the sample that is CPLEN samples 
before the start point determined in step 408. In step 414, the 
coefficients are returned for use in the TEQ filter. 

While, as an alternative, it is easy to compute all possible 
sets of TEQ coefficients and then choose the set that pro- 
vides the best performance, such a technique is time con- 
suming. On the other hand, the technique described herein 
for determining the TEQ coefficients is much less burden- 
some. 

The above- disclosed algorithm may be implemented in 
any of the known ways of implementing TEQ filters, 
including, but not limited to, a programmed microprocessor, 
a micro controller, a state machine, combinatorial logic or 
any combination of the aforementioned. Thus, in accordance 
with the present invention, an efficient method and apparatus 
for calculating the coefficients for a time domain equaliza- 
tion filter is provided along with an efficient frequency 
domain apparatus and method for symbol synchronization. 

Having thus described a few particular embodiments of 
the invention, various alterations, modifications, and 
improvements will readily occur to those skilled in the art. 
Such alterations, modifications and improvements as are 
made obvious by this disclosure are intended to be part of 
this description though not expressly stated herein, and are 
intended to be within the spirit and scope of the invention. 
Accordingly, the foregoing description is by way of example 
only, and not limiting. The invention is limited only as 
defined in the following claims and equivalents thereto. 

I claim: 

1. Amethod of generating coefficients for performing time 
domain equalization on a digital multitone signal compris- 
ing a plurality of symbols received from a transmitter over 
a communication channel, said method comprising the steps 
of: 

(1) synchronizing said received digital multitone signal to 
a corresponding signal that was transmitted from said 
transmitter prior to inclusion of any channel response; 
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where V^V^R^" 1 , 

V r is a lower triangular matrix whose main diagonal 

elements are all l J s, and 
D^ v is a diagonal matrix, 

3. The method as set forth in claim 2 wherein step (2) 
comprises the steps of: 

(2.1) calculating each of V r and D tm , once; 

(2.2) storing V r and D inv in a memory; and 

(2.3) retrieving V T and when needed for determining 

4. The method as set forth in claim 2 wherein step (2) 
further comprises: 

(2.4) calculating VR Ay once and using that calculated 
value for each occurrence of VR A> , in tbe equation. 

5. The method as set forth in claim 2 wherein step (2) 
comprises: 

(2.5) solving the equation for a. opt by employing a 
Cholesky decomposition technique. 

6. The method as set forth in claim 1 wherein said method 
is used in connection with a digital multitone receive archi- 
tecture that utilizes a cyclic prefix of length L and wherein 
step (2) provides a coefficient wherein M^L, where M is the 
impulse response length of a shortening channel utilized in 
said method. 

7. The method as set forth in claim 6 wherein L-M. 

8. The method as set forth in claim 7 wherein step (2) 
comprises assuming that an impulse response of a channel of 



10 



(2) determining time domain equalization coefficients, 
SL opt , for a TEQ filter via 

r 

where 

iU*)^M«M»+*)}^(-*) 
^*>^««M»+*)}^(-*) 
V*>^M»M»+*)}^(-*) 

E{ } in the equations above represents a generalized 

averaging function or an ensemble average 

and 



said received digital multitone signal has a transfer function 
of: 



r = E{y{n)W{> 



W(nHx(n) x(n-l) x(n-2) . . . x(n-(I^l))y(n-l) y(n- 

2)...y(n-(K-l)f 
L is the length of a cyclic prefix 
K is the length of the TEQ filter 
T represents a transpose function; and 
(3) time domain equalizing said received signal using said 

coefficients determined in step (2). 
2. The method as set forth in claim 1 wherein step (2) 
comprises solving for a opt using one of 



25 



9. The method as set forth in claim 6 wherein said time 
10 domain equalization is performed in connection with an 

asynchronous digital subscriber line. 

10. The method as set forth in claim 9 wherein said time 
domain equalization is used in connection with a commu- 
nications system in accordance with ITU-T G.992.2 speci- 

15 fications. 

11. The method of claim 1 wherein step (1) comprises: 

(1.1) cross correlating the received signal with the trans- 
mitted signal to generate an impulse response of the 

20 channel; 

(1.2) determining a maximal energy window of said 
impulse response and a range of said channel response; 

(1.3) determining a noise floor of said channel impulse 
response; and 

(1.4) determining a start point of said channel impulse 
response using said noise floor. 

12. The method as set forth in claim 11 wherein step (1.2) 
comprises passing said impulse response through an energy 

30 detector with a window of size L samples. 

13. The method as set forth in claim 12 wherein step (2) 
is performed starting with a sample of said received signal 
that is L samples before said start point. 

35 . 14. The method as set forth in claim 11 further comprising 
the steps of: 

(1.5) converting said received signal to frequency domain 
before step (1-1); and 

(1.6) converting said channel impulse response from 
frequency domain to time domain after step (1.1). 

15. A method of recovering a signal that was received 
after, being transmitted in digital multitone with a cyclic 
prefix, said method comprising the steps of: 
45 (1) time domain equalizing said signal; 

(2) removing said cyclic prefix from said signal; 

(3) frequency domain equalizing said signal; 
wherein step (1) comprises convolving said signal with 

50 l+a(h) wherein a(n) is generated via 

or 

55 

a^RjR^R^-R^T^-RjR^rJ-, 

where 

60 «^)»E{x(n)rM}^(-i) 
R x £k)=E{x(n)y(n + k)}=R y X-k) 
RyW-EiyWyin+ty-R^-k) 



65 



E{ } in the equations above represents a generalized 
averaging function or an ensemble average 
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and 



--E{y{nmn)} = 



where 

W(n)=[x(n) x(n-l) x(n-2) . . . x(n-(L-l))y(n-l) 

y(n-2) . . . y(n-(K-l) ] T 
L is the length of a cyclic prefix 
K is the length of a TEQ filter 
T represents a transpose function. 
16. The method as set forth in claim 15 wherein step (1) 
comprises solving for a apt using one of 



10 



12 



(4.5) converting said signal to frequency domain before 
step (4.1); and 

(4.6) converting said channel impulse response from 
frequency domain to time domain after step (4.1). 

28. A receiver for a digital multitone modem comprising: 
means for synchronizing a received digital multitone 
signal to a corresponding signal that was transmitted 
from a transmitter prior to inclusion of any channel 
response; 

means for determining, for said received signal, time 



domain equalization coefficients, a, 
via 



for a TEQ filter 



(0«^) r i) to VH^-il>^r I -((W^,) , 3) i((|( W 1 )-, 

where V^y-R^" 1 , 

V T is a lower triangular matrix whose main diagonal 

elements are all l's, and 
D /nv is a diagonal matrix. 

17. The method as set forth in claim 16 wherein step (1) 
further comprises: 

solving one of said equations for & opt by employing a 
Cholesky decomposition technique. 

18. The method as set forth in claim 15 wherein said 
cyclic prefix is of a length, L, and wherein step (1) shortens 
a channel of said signal by M, where M^L. 

19. The method as set forth in claim 18 wherein L=M. 

20. The method as set forth in claim 19 wherein step (1) 
further comprises assuming that an impulse response of a 
channel of said signal has a transfer function of: 



: Mz) " 



20 



25 



35 



21. The method as set forth in claim 18 wherein said 
signal is transmitted via an asynchronous digital subscriber 
line. 

22. The method as set forth in claim 21 wherein said 
method is used in connection with a communications system 
in accordance with ITU-T specification G. 992.2. 

23. The method as set forth in claim 15 further comprising 
the steps of: 

(4) determining a starting point for calculating a. opr 

24. The method of claim 23 wherein step (4) comprises: 

(4.1) cross correlating the received signal with the trans- 
mitted signal to generate an impulse response of the 
channel; 

(4.2) determining a maximal energy window of said 
impulse response and a range of said channel response; 

(4.3) determining a noise floor of said channel impulse 
response; and 

(4.4) determining a start point of said channel impulse 
response using said noise floor. 

25. The method as set forth in claim 24 wherein step (4.2) 
comprises passing said impulse response through an energy 
detector with a window of size L samples. 

26. The method as set forth in claim 25 wherein step (1) 
is performed starting with a sample of said received signal 
that is L samples before said start point. 

27. The method as set forth in claim 24 further comprising 
the steps of: 



60 



where 

R xy {k)=E{x(n)y{n+k)}^ yx (-k) 

E{ } in the equations above represents a generalized 

averaging function or an ensemble average 

and 



r = E{y(n)W(n)) 



where 



=11 



45 



50 



W(n)-[x(n) x(n-l) x(n-2) . . , x(n-(L-l))y(n-l) y(n- 

2) . . . y(n-(K-l)f 
L is the length of a cyclic prefix 
K is the length of the TEQ filter 
T represents a transpose function; and 
means for time domain equalizing said received signal 

using said coefficients. 
29. The receiver as set forth in claim 28 wherein said 
means for determining time domain equalization coefficients 
comprises solving for z. opt using one of 



55 



where V^V-R^ 1 , 

V r is a lower triangular matrix whose main diagonal 
elements are all l*s, and 
is a diagonal matrix. 

30. The receiver as set forth in claim 29 wherein said 
means for determining time domain equalization coefficients 
solves the equation for a^, by employing a Cholesky 

65 decomposition technique. 

31. The receiver as set forth in claim 28 wherein said 
receiver is used in connection with a digital multitone 
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15 



34. The receiver as set forth in claim 31 wherein said 
receiver is used in connection with an asynchronous digital 
subscriber line. 

35. The receiver as set forth in claim 34 wherein said 20 
receiver is used in connection with a communications sys- 
tem in accordance with ITU-T specification G. 992.2. 

36. The receiver as set forth claim 28 wherein said means 
for synchronizing comprises: 



14 



receive architecture that utilizes a cyclic prefix of length L 
and wherein said means for determining time domain equal- 
ization coefficients provides a coefficient wherein M^L, 
where M is the impulse response length of a shortening 
channel utilized in said method. s 

32. The receiver as set forth in claim 31 wherein L-M. 

33. The receiver as set forth in claim 32 wherein said 
means for determining time domain equalization coefficients 
assumes that an impulse response of a channel of said 10 
received digital multitone signal has a transfer function of: 



means for cross correlating the received signal with the 
transmitted signal to generate an impulse response of 
the channel; 

means for determining a maximal energy window of said 
impulse response and a range of said channel response; 

means for determining a noise floor of said channel 
impulse response; and 

means for determining a start point of said channel 
impulse response using said noise floor. 

37. The receiver as set forth in claim 36 wherein said 
means for determining a maximal energy window comprises 
means for passing said impulse response through an energy 
detector with a window of size L samples. 

38. The receiver as set forth in claim 36 wherein said 
means for synchronizing comprises: 

means for converting said received signal to frequency 
domain before cross correlating the received signal 
with the transmitted signal; and 

means for converting said channel impulse response from 
frequency domain to time domain after cross correlat- 
ing the received signal with the transmitted signal. 
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